Modulated currents in open nanostructures 
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We investigate theoretically the transport properties of a mesoscopic system driven by a sequence 
of rectangular pulses applied at the contact to the input (left) lead. The characteristics of the 
current which would be measured in the output (right) lead are discussed in relation with the 
spectral properties of the sample. The time-dependent currents are calculated via a generalized 
non-Markovian master equation scheme. We study the transient response of a quantum dot and of 
a narrow quantum wire. We show that the output response depends not only on the lead-sample 
coupling and on the length of the pulse but also on the states that propagate the input signal. We 
find that by increasing the bias window the new states available for transport induce additional 
structure in the relaxation current due to different dynamical tunneling processes. The delay of the 
output signal with respect to the input current in the case of the narrow quantum wire is associated 
to the transient time through the wire. 
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I. INTRODUCTION 

Time-dependent transport measurements at nanoscale 
provide important insight into the intrinsic properties of 
semiconductor structures like relaxation and dephasing 
timesi and play a crucial role in single-shot spin read-out 
schemesi^ Consequently transient response of quantum 
structures to pulsed signals gains interest from experi- 
mental point of view. Recently Naser et al^ measured 
the time-dependent current through a quantum point 
contact when a pulse generator is coupled to the input 
lead. The output signal was measured for different am- 
plitudes and rise times of the pulse. The decay of the out- 
put current was roughly exponential and it was suggested 
that such a device could be used as a microwave circuit 
element. In another recent experiment Lai et alA per- 
formed transient current measurements for a Ge quantum 
dot when trapezoidal voltage pulses were applied to one 
electrode. The transient current mimics the pulse shape 
and the extracted time-dependent occupation number in 
the system follows an accumulation/depletion cycle. 

From the theoretical point of view the transient cur- 
rents have been calculated by various methods. Ste- 
fanucci et al^^ combined the Keldysh formalism and 
the density functional theory to study the response of 
a nanostructure to a time-dependent bias between the 
leads. In the partitioning approach the transient currents 
induced through a few-level quantum dot by a modula- 
tion of the contacts between the leads and the sample was 
calculated in Ref0 and a scattering theory of the time- 
dependent magnetotransport in a long quantum wire was 
proposed in RefH. 

These studies were essentially focused on understand- 
ing the transient regime and the onset of the steady-state. 
Therefore the time-dependent driving used in the calcula- 
tions did not describe pulsed signals, but rather an initial 
switching stage followed by a constant, time- independent 
value. In particular sudden coupling or smooth coupling 



lead to qualitatively different output signals. Another 
problem considered in the time-dependent transport cal- 
culations is the quantum pumping i^ d°i^^d^ In this case 
one obtains averaged currents through an unbiased sys- 
tem which is perturbed by two time-dependent potentials 
oscillating out-of-phase. 

In this work we discuss transport calculations for a 
mesoscopic sample driven out of equilibrium by both a 
constant bias applied between the leads and a fast os- 
cillating signal applied at the contact between the lead 
and the sample. If the signal has a rectangular shape the 
setup is in some sense similar to the pump-and-probe 
configuration used in the transient spectroscopy experi- 
ments of Fujisawa et alX The time-dependent signal was 
applied on the sample and the main aim was to extract 
the spin relaxation time by pushing one excited state into 
the transport window during the pulse. 

Here we discuss the transient response of the system 
from a different angle: If a sequence of rectangular pulses 
that modulate the coupling to the left lead is viewed as an 
input signal then one can study the propagation of this 
signal trough the system and the corresponding current 
in the right (i.e. output) lead. Our problem is therefore 
closer to the experiments by Naser et al^ and Lai et alA 
Besides the very ambitious goal of assembling quantum 
dot structures in complex mesoscopic circuits operating 
like quantum gates we believe there is another important 
motivation for such a study. The transient response of 
the sample to the modulation of the contact is a conse- 
quence of the internal electron dynamics which depends 
crucially on the electronic states participating in trans- 
port. The point is that by changing the bias applied 
on the system one selects different states in the trans- 
port window and then the output current may carry im- 
portant informations about the electron dynamics in the 
sample. Even if the level structure of the sample may be 
known from other type of measurements, the tunneling 
rates associated to each quantum state or the propaga- 
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tion properties of the 'orbitals' are not easily understood. tion): 

We aim at describing the following transport exper- 
iment through a mesoscopic structure: the contact be- 
tween the sample and the left lead opens and closes pe- 
riodically by applying rectangular pulses on the metallic 
gates that define the contact region. At the same time the 
contact between the sample and the right lead gradually 
opens. Then each time the left contact closes the elec- 
trons in the sample can only escape into the right lead. 
These relaxation processes of various states in the sam- 
ple depend mainly on the corresponding tunneling rates 
which in turn are given by the coupling of the states to 
the output contact. In real samples the states that par- 
ticipate in transport have different tunneling coefficients 
and therefore it is not obvious that the relaxation cur- 
rent follows a simple exponential decay. We show that 
the transient response of the sample is in general more 
complex. Sometimes the output current may look expo- 
nential, but often may also carry a fine structure reflect- 
ing the presence of several quantum states propagating 
through the sample and various relaxation processes cor- 
responding to transitions between these states and the 
right lead. 



A recent attempt to describe rectangular pulses in 
quantum dots was made by Oh et al^ Their sample 
model was a two-level quantum dot and for transport cal- 
culations they used rate equations with time-dependent 
tunneling coefficients. Instead, our calculations are done 
using the generalized non-Markovian master equation 
(GME) approach presented in one of our recent worksJ^ 
We take into account the geometry of the system and its 
spectral properties which determine the tunneling coeffi- 
cients and therefore the currents driven by the external 
bias. The application of the GME method to quantum 
transport received attention in the last years in the con- 
text of transient currents, measurement theory or full- 
counting statistics or coherent control of transport (see 

e.g. Refs.lillilliKaiilliiliaiil). 

The content of the paper is organized as follows: the 
formalism is briefly explained in Section II, the numeri- 
cal simulations are reported in Section III, and the main 
conclusions are included in Section IV. 



II. THE GENERALIZED MASTER EQUATION 
METHOD 



In this Section we introduce the Hamiltonian of the 
system, the underlying notations, and we summarize the 
main equations derived in our previous work Reffisl. We 
consider a mesoscopic system which is coupled to two 
leads at t = and described by the following second- 
quantized Hamiltonian (/i.e. denotes Hermitian conjuga- 
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l=L,R'' n 
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where c|^;, Cqi, and dj^, d„ are creation and destruction op- 
erators for electrons with momentum q and energy (q) 
in the lead /, and with energy En in the sample, respec- 
tively. The labels L and R denote the left and right lead. 
The third term is the tunneling Hamiltonian and con- 
tains the time-dependent switching functions x'(^) ^.nd 
the coupling matrix elements T^^^ associated to each 

pair of states {ipq ^ ^ , (j)n\ from the leads and the sam- 
ple. The pulses applied at the contact region between 
the left lead and the sample are simulated by the func- 
tion XL which by construction has a rectangular shape. 
Although our formalism could also be implemented for 
a continuous model (see Ref.[23l) here we use a lattice 
model for which the matrix elements are given by (see 
Ref. dll)): 



(2) 



where Vi is the coupling strength between the lead I and 
the sample, 0; is the site of the lead I which couples to 
the contact site ii in the sample. The eigenfunctions of 
the sample and the corresponding energies are 
numerically computed while those in the leads, ?/;^ and 
e'(g) are known analytically: 
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\q)^2tLC0sq, (3) 



being the hopping constant on leads. We emphasize 
that even in this simple lattice model the coupling ma- 
trix elements introduced in Eq.([2]) depend both on the 
energy and of the localization of the sample states. A 
similar model for the transfer Hamiltonian was proposed 
by Maddox et alM 

The statistical operator of the open quantum sys- 
tem is denoted by W and solves the Liouville equa- 
tion: iW{t) = [H{t),W{t)], with the initial condition 
W{t = 0) = plPrPs- This means that before the cou- 
pling the sample is described by the statistical operator 
ps defined just below, and the leads are characterized 
by equilibrium distributions pl,r with different chemical 
potentials p,L > Pr- 

The many-body states of the system arc described 
by the sequence of occupation numbers of the single- 
particle states {4>n} for the isolated system. We shall 
denote the many-body states by Greek letters, i.e. {ly) = 
^2, .., and by the occupation number of the 
n-th single particle state. If the initial state of the dis- 
connected sample is \vo) then = |i'o)('^o|- For example, 
if two electrons are situated on the lowest levels at t < 
we have Ivq) = |110....). 
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Following the main lines of the superoperator 
method^^i^ we take the partial trace of W over the Fock 
space of the leads and end up with a master equation for 
the reduced density operator p{t) = TTLTrfiW{t) with 
the initial condition p(t = 0) = ps up to second order in 
the tunneling Hamiltonian: 




where we have introduced the operators (see Ref. for 
details): 

J to 

^i/3('7)-E^n«H4l/3)- (5) 

n 

It is clear that T^p{q) describes the 'absorption' of elec- 
trons from the leads to the system and changes the many- 
body states of the latter from /? to a. Observe that 
T2f^{q) 7^ only if the number of electrons in the many- 
body states a and /3 differ by one. /; denotes the Fermi 
function in the lead I. The difference between the chem- 
ical potentials defines the bias applied across the sample 
eV = PL — PR- Observe also the presence of loss and 
gain terms in Hqi . 

We also define a set of relevant states located in the 
energy window [£^min , £^max] , where -Emin < pr < pl < 
Emax- This 'active' window includes only those states 
of the sample that are relevant to the transport. More 
precisely, i^min is chosen such that the levels with lower 
energy are fully occupied both prior to the couphng of the 
leads, i.e. for t < 0, and also after the coupling began, 
at < > 0, in the presence of the bias. Similarly, -Bmax is 
selected such that all states with higher energy are per- 
manently empty. Consequently the states outside this 
energy window do not contribute to the current. Based 
on this picture we conclude that it is sufficient to com- 
pute an 'effective' reduced density matrix by taking into 
account only those many-body configurations resulting 
from the single particle states within the active window. 
Also for the simplicity of notation we shall specify in 
the many-body states only the occupation numbers of 
the single-particle states within the active window. Of 
course, the validity of this truncation should be checked 
in the numerical simulations by gradually enlarging the 
active window [£^minj £'max] until the calculated currents 
become stable. 

The time evolution of the charge residing in the active 
region is related to the diagonal elements of the reduced 



density matrix: 

{Qs{t))=Y.T.'n{'^\pit)W)- 

n V 

Using the GME, Eq. ([5]), one can easily identify the con- 
tribution of each level n to the currents in the left and 
right lead: 

■h = E Jlu 

n 

•^'^" = -^E*« / dqx\t){u\[Ti,n,i{t)] + h.c.W) 

(6) 

where i^j = 0,1 is the occupation number of the n-th 
single particle state inside the active window. The sign 
of the net currents in the leads is positive if it is oriented 
from the left to the right, i.e. > if the electrons flow 
from the left lead towards the sample and J^j > if they 
flow from the sample towards the right lead. During the 
transient regime the sign of the net currents may change 
in time. 

The GME is numerically solved through the Crank- 
Nicolson method (see the details in Ref. (|l3l)). Through- 
out this work the Coulomb interaction effects are not 
considered; further discussion on this point is given at 
the end of Section III. 



III. NUMERICAL RESULTS 

The first system we consider is a rectangular two- 
dimensional lattice having 25 x 20 sites. The correspond- 
ing Hamiltonian has 500 eigenvalues £"„ and eigenvectors 
(^„, n = 1, 2, 500, as many as the number of sites. The 
energy unit is given by the hopping parameter in the 
sample t_D = fi,'^/2m*a^, where a is the lattice constant 
and TO* is the electron effective mass in GaAs. The spec- 
trum of the isolated sample is contained in the range 
(— 2t£), 2t£)). For a = 8 nm our lattice describes a 200 
nm xl60 nm sample. We fix kT = 0.1 meV which cor- 
responds to a low temperature T = 1.15 K. The sites of 
the system are denoted by i and are specified by the pair 
of coordinates (ix,iy)- 

In order to describe the gradual couphng of the leads to 
the sample and the periodic modulation of the left con- 
tact we use specific coupling functions Xi{t)- The cou- 
pling to the leads begins at t = and evolves in time 
like f{t) = 1 — 2/{e''* + 1). The parameter 7 defines the 
smoothness of the coupling. After some time to > the 
coupling to the left lead is turned off and on periodically. 
These pulses are analytically defined by combining the 
functions f{t) and 1 — f{t). We denote the pulse length 
by Tp. We think this kind of time-dependent perturba- 
tion is a reasonable model of an experiment in which 
the sample is smoothly coupled to the leads with differ- 
ent chemical potentials, possibly reaching an equilibrium 
state before the pulses begin to act. 
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FIG. 2: (Color online) (a)-(c)Total currents in the input (left) 
and output (right) leads for the 25 x 20 sites system for differ- 
ent pulse lengths: (a) Tp — 7.5 ps, (b) Tp = 15 ps, (c) Tp = 30 
ps. The bias window covers only two single particle states. 
In Figs, (b) and (c) we show the currents for longer times 
than in figure (a) in order to capture the same number of 
pulses. The signal applied on the input contact is also qual- 
itatively shown with a dashed line. The coupling strength 
Vl = Vr= 4, = -3.65 and = -3.7. 



FIG. 1: (Color online) The localization probabilities asso- 
ciated to the four states included in the active window. 
£10 = -3.68 (fc = 1), Ell = -3.67 {k = 2), E12 = -3.63 
{k = 3) and £13 = -3.62 {k = 4) (in to units). The left lead 
is located near the lower left corner of the sample and the 
right lead coupled close to the upper right corner. 



The currents depend strongly on the placement of the 
contacts. In this example the two leads are attached 
at diagonally opposite corners of the sample. We first 
choose the chemical potentials /i^ = —3.65 and /i/j = 
—3.60. Then our sample has nine energy levels below /i/j, 



and as we checked they do not contribute significantly to 
the transport. The chosen active region contains two 
states within the bias window which are £^10 and En 
and two more states above fiL, which are E12 and E13. 
Instead of the label n ~ 10, 11, 12, 13 for the states in the 
active window (in the present example), we will also use 
the label fc = 1, 2, 3, 4 respectively. 

In Figs. l(a)-(d) we show the site occupation probabil- 
ities |(/)„(i)p for the active states. We emphasize that the 
middle state with n = 12 or k = 3, Fig. 1(c). is weakly 
coupled to the leads while the other ones have a larger, 
but still moderate coupling. 
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FIG. 3: (Color online) The total charge accumulated on the 
two states from the bias window for the three pulse lengths 
considered in Figs. 2(a)-(c). 



Figs. 2(a)-(c) show the time-dependent total currents 
in both leads for different pulse lengths. We also indi- 
cate qualitatively the pulsed signal as given by the func- 
tion XLit) (arbitrary units are used on the corresponding 
axis) . In the beginning the system is coupled to the leads 
and both contacts are kept open for a time to. Then the 
the left contact is modulated by pulses while the other 
one is left open. As the sample dynamics adapts to the 
input signal a periodic regime is already established after 
two pulses. The first observation is that for very short 
pulses (Tp = 7.5 ps in Fig. 2(a)) the output current 
has a triangular shape although the modulating signal is 
rectangular. In addition Jji does not vanish when the left 
lead is disconnected, resembling the charging-rclaxation 
characteristic of a capacitor. The charge is first absorbed 
from the left lead when the contact is switched on, and 
then only partially expelled into the right lead, when the 
contact to the left lead is off, i.e. even in the absence of a 
driving bias. The decay of the output signal looks expo- 
nential, but after one complete cycle, i.e. before the left 
contact opens again, the magnitude of the output cur- 
rent in the right lead is still considerable. The current in 
the left lead Jl decays much faster when the contact is 
turned off, following closer the modulation potential at 
the contacts. We notice that right after the left contact 
opens the current is injected in the sample quite fast but 
the current in the right lead increases slower. 

When the pulse length increases the shape and the 
amplitude of the output current change considerably. Jn 
reaches maxima even before the pulse is turned off and 
remains almost constant during the second half of the 
pulse, while at the same time the input current decrease. 
However, at Tp ~ 15ps the output current still shows 
an exponential decay and does not reproduce the input 
signal. A different response to the pulse train is obtained 
at Tp = 30ps. In this case the output current decays first 
exponentially and then stays flat until the left contact 
opens again. 
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FIG. 4: (Color online) (a)-(c)Total currents in the input (left) 
and output (right) leads for the 25 x 20 sites system for differ- 
ent pulse lengths: (a) Tp = 7.5 ps, (b) Tp = 15 ps, (c) Tp = 30 
ps. The bias window covers four single particle states. The 
signal applied on the input contact is also given in arbitrary 
units- dashed line. The coupling strength Vl = Vr = 4. 
yLh = —3.6 and ^lr — —3.7. 



Comparing the behavior of the input currents in all 
three cases we see that in the first half of the 30 ps pulse 
the current is similar to the one that develops in the full 
length of the 15 ps pulse. The saw-tooth profile of the 
current in Fig. 2(a) is also present in the first half of the 
15 ps pulse. This is due to the fact that as long as the 
sample is coupled to both leads it has the transient be- 
havior which is already observed during the initial charg- 
ing time. It is also clear that if the pulse length is too 
short the left lead does not feed enough charge to the 
sample in order to maintain a constant output current. 

An estimate of the pulse length which generates an 
output current with almost a rectangular shape can be 
taken from the first transient period (t < to); that pulse 
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FIG. 5: (Color online) Partial currents in the output (a) and 
input (b) leads for the 25 x 20 sites system for Tp — 7.5 ps 
and = —3.6 and hr = —3.7. In this case the bias window 
covers four states, (c) The output partial currents for the two- 
level configuration /iz, — —3.65 and — —3.7. The coupling 
strength is Vl = Vr = 4. 



length should be at least equal to the time at which the 
output current become nearly equal. 

In Fig. 3 we compare the total charge accumated on 
the two states within the bias window for the three pulse 
lengths. As the pulse length increases more charge is 
transferred through the system and therefore the output 
current increase. One should notice that for the 30 ps 
pulse the charge first relaxes exponentially but then al- 
most linearly. Since the current is essentially the deriva- 
tive of the charge with respect to time this means the 
current in the right lead is first exponential and then 
constant, as we have already learned from Fig. 2(c). 

In Figs. 4(a)-(c) we show the currents obtained from 
the same sample when the chemical potential of the 
left lead is pushed to fiL — —3.6 and thus two more 
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FIG. 6: (Color online) (a) The occupation numbers of states 
from the bias window for a pulse length Tp — 7.5 ps, fJ.L = 
—3.6 and hr = —3.7. (b) The same in the two level case 
= —3.65 and /li? = —3.7. 



states enter into the bias window, i.e. those with fc = 3 
and and fc = 4. When comparing to Fig. 2 we observe 
additional 'shoulders' developing in the output current. 
These shoulders are produced by the new states included 
in the bias window during the relaxation into the right 
lead. The relaxation processes depend on the coupling 
between the states and the lead. Since one of the two 
new states is only poorly coupled to the leads, i.e. the 
one withe fc = 3 (see Fig. 1(c)), the shoulder is actually 
produced by the new state with fc = 4. 

In the following we compare the partial currents in 
the right lead for the two states and four states config- 
urations in the case of the ultrashort 7.5 ps pulse. The 
partial currents Ji^k are calculated with Eq.(IB]). Fig. 5(a) 
gives the output currents and confirms that the 3rd state 
does not contribute significantly to the transport. Both 
Jr^i and Jfi^2 increase during the entire pulse length, al- 
though with a lower slope in the second part. Rather 
surprisingly, in this time interval Jr^4 decreases. In the 
relaxation interval the situation is the opposite: Jf{,i and 
Jr.,2 are monotonously decreasing while the current of 
the 4th state considerably increases in the second half of 
the relaxation interval. The currents entering the sample 
from the left lead (see Fig. 5(b)) also have interesting fea- 
tures: Jl^i and Jl,2 rise suddenly to a maximum value, 
while Jl^4 increases slower and does not reach a maxi- 
mum within the pulse. This suggests that the two lowest 
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states absorb quickly more charge from the left reservoir. 
By looking at the occupation numbers shown in Fig. 6(a) 
one convinces himself that this is indeed the case. The 
occupation of the 4th level increases much slower than 
iVi and N2. 

The behavior of Ji^4 and Jr^^ can be explained by the 
dynamic tunneling processes described by the gain and 
loss terms of the solution of the master equation via the 
matrix 11^;. Suppose, for example, that at time s < < an 
electron tunnels from the input lead into the lowest state 
of the bias window, fc = 1 . At instant t the same electron 
may tunnel out into the right lead. But it is also possible 
that this electron remains in the sample, while another 
electron, from another state, say the 4th, tunnels back 
into the left lead; escaping into the left lead from the 4th 
state is more likely than in the right lead because the 
chemical potential is much closer to Ei^ {k = 4) than 
/ifl. Repeated tunneling from highest state of the sample 
towards different leads and back into the sample lead 
to the charging of the lower states at the expense of the 
higher ones. Consequently the output currents associated 
to the lower states 1 and 2 are increasing, see Fig. 5(a), 
whereas the output current J_r4 decreases. Note also 
that the net input current Jl,a is smaller than J^^i and 
Jl,2 because there are more electrons tunneling out from 
the level k = A. 

In contrast, when the pulse is turned off, the charge on 
the 4th level leaves the sample only via the output lead 
and therefore Jn^ increase as shown in Fig. 5(a). On the 
other hand the tunneling processes from the left lead to 
the lowest levels are switched off and so Jfl i and J/f^2 
decrease on the relaxation interval. 

For comparison we also show in Fig. 5(c) and Fig. 6(b) 
the output currents for the two-level configuration (fc = 
1,2) and the occupation numbers of the two levels. For 
the two-level configuration the partial occupation num- 
bers display similar charging/relaxation shapes. This 
happens because these two particular states are equally 
coupled to the leads and hence the input signal propa- 
gates almost identically through both states towards the 
right lead. Remarkably, a similar behavior of the time- 
dependent occupation number was obtained from the ex- 
perimental data by Lai et alA 

In order to get more informations on the relevant tun- 
neling processes we have analyzed the diagonal elements 
(i.e. the populations) of the reduced density operator. 
It useful to introduce a shorter notation for the many- 
body states by interpreting the occupation numbers as 
decimal numbers written as binary strings (but reading 
them in the reverse order, and adding 1). For example 
|1) = 10000), |2) |1000), |3) = 10100) etc. 

Fig. 7(a) shows the populations corresponding to the 
single-particle sector of the Fock space, that is, the prob- 
abilities for all the configurations containing only one 
electron. We also give in Fig. 7(b) the most relevant two- 
particle configurations. The vacuum state |0000) has the 
largest probability pn which is not shown. At t = it is 
1 because the sample is initially empty, then during the 
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FIG. 7: (Color online) The relevant diagonal elements of the 
reduced density operator for the 25 x 20 sites system cal- 
culated with four single particle states in the active region, 
(a) The diagonal elements of the reduced density operator 
corresponding to single-particle configurations |2) = |1000), 
|3) = lOlOO), |5) = lOOlO), |9) = lOOOl). (b) The most rel- 
evant two-particle configurations |4) — |1100), |6) — |1010), 
|7) = |0110> |10) = |1001). HL = -3.6 and hl = -3.7. 



charging period it drops to about 0.6, but it increases 
back to about 0.7 during the relaxation interval. From 
Fig. 7(a) we infer that i) the configurations |1000) and 
|0100) are the most probable in the transient regime; ii) 
the corresponding probabilities P22 and ^33 decay expo- 
nentially (qualitatively speaking) during the relaxation, 
but pgg which is the probability of the state |0001) does 
not decay exponentially. This state is even stable for 
some time in the relaxation regime, while the probabil- 
ities for the states |1000) and |0100) decrease. This is 
another way of seeing that the 4th state relaxes later 
into the right lead, iii) P55 is only slowly increasing and 
has only small oscillations due to the time-dependent sig- 
nal, due to the weak coupling to the leads. The two- 
particle configurations shown in Fig. 7(b) have very all 
small probability because the system is poorly charged 
during the ultra-short pulse considered here. 

For completeness we give in Figs. 8(a)-(c) the off- 
diagonal elements of the reduced density matrix, known 
as " coherences" , which are related to the transitions be- 
tween different states. For example P29 (Fig. 8(a)) de- 
scribes the process in which electrons can tunnel from 
the 4th state into the leads and then back to the sam- 
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FIG. 8: (Color online) (a)-(c) The imaginary and real parts 
of the off-diagonal elements of the reduced density matrix 
which couple various states (see the discussion in the text). 
Note the multiple oscillations in any half-period and the small 
amplitude for p29. The dashed (blue) line represents the signal 
applied on the input contact. Tp = 15 ps, /xl = —3.6 and 
^.L = —3.7. 



pic on the 1st state. We see that the real and imaginary 
parts of the coherences have oscillations and change sign. 
Moreover, the oscillations are more complex if the matrix 
element implies single particle states separated by other 
energy levels. Indeed, by comparing Rep23 and Re/929 we 
see that the latter displays two minima and two max- 
ima in the charging time interval, while Rep23 has only 
one minimum. This feature suggests again that electrons 
make transitions via intermediate single-particle states 
as follows: initial state intermediate state — s- leads 
final state. We have checked that coherences between 
states with only one intermediate state in between (e.g. 
Re/939) have only one maximum and one minimum. 
Another remark is that the oscillations of p23 have a 



FIG. 9: (Color online) (a) The total input and output currents 
for a 5 X 20 sites system with three states within the bias 
window. Jij is delayed w.r.t. to Jl- The pulse is rather 
squared (the dashed (green) line), with the rise time defined 
by the parameter 7=1. (b) The occupation numbers of the 
states within the bias window. Note that the lowest states 
relaxes at a slower rate due to a smaller coupling to the leads, 
(c) The total input and output currents for a smooth pulse 
(7 = 0.25). Tp = 15 ps,/iL — —3 and = —3.6. 



higher amplitude than the other two coherences, which 
is consistent to the fact that the two configurations 
|2) = |1000) and |3) = |0100) are the most probable ones 
in the stationary regime. Fig. 8(c) shows that the 3rd 
level is very weakly coupled to the level below, which is 
expected since that state is poorly coupled to the leads. 
By increasing the pulse width the number of oscillations 
during the pulse generally increases as more and more 
transitions may occur. Also, two-particle configurations 
will develop, while the probability of the single-particle 
states decrease (not shown). It should be mentioned here 
that in the non-Markovian approach the coherences also 
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contribute to the currents because they are coupled to 
the diagonal elements of p. 

We continue with similar results for a different sample, 
a narrow 40 nm x 160 nm quantum wire described by a 
5 X 20 = 100 sites lattice. The system has now 100 eigen- 
states and it is connected to two three-channel leads on 
both narrow sides. For simplicity we consider uncoupled 
channels, which means the electrons cannot jump be- 
tween the channels. The chemical potentials of the leads 
are = ~3 and ^.j^ = —3.6; they cover three states now, 
^3 = -3.53, Ei = -3.38 and = -3.19 (or /c = 1, 2, 3 
respectively) which we also consider as the active win- 
dow. The coupling of these states to the leads is much 
stronger than in the previous case. Fig. 9(a) shows the 
input/output characteristics for a pulse of length Tp ~ 7.5 
ps. We see that the total currents become equal in the 
second half of the pulse. This means that the station- 
ary state is practically reached before the left contact 
is switched off. After that the output current is delayed 
with respect to the input current. This delay corresponds 
to the propagation of electrons along the sample. We 
have checked this by estimating the time needed by one 
electron having the energy E4 = 2tL cosg to travel along 
the wire; we obtain that the traveling time is around 2 
ps which agrees with the numerical results. 

The input current has a sharp transient maximum as 
the left contact opens, but this behavior is not trans- 
ferred to the output current which is rather flat. Again 
the relaxation input current has a 'step' after a fast de- 
cay. The occupation numbers of the three states within 
the bias window are shown in Fig. 9(b) and explain the 
origin of the step. The two higher states depopulate quite 
fast in the first half of the relaxation interval, but at some 
point their depletion suddenly slows down. This suggests 
that the relaxation rates are governed by time-dependent 
tunneling rates depending on the coupling to the leads, 
but also on the occupation number of the states. Re- 
mark that the bias window becomes empty before an- 
other pulse rises up. 

The rise time of the pulse can be varied in experiments 
(see Ref. i). In our model the rise time depends on the 
smoothness parameter 7 used to create the pulses. In 
Fig. 9(c) we show the total input and output currents for 
a pulse with a larger rise time corresponding to 7 = 0.25. 
The transient maxima in the input current softens and 
so does the additional shoulder in the relaxation current. 
Note however that during the pulse the output current 
takes the same value (around 60 pA) and it is still delayed 
w.r.t. to the input current. 

Finally we would like to comment on the Coulomb in- 
teraction which is neglected here. The numerical results 
show that the occupation number of the levels within the 
active window is quite low for the rather large samples 
considered (see e.g. Figs. 6(a) and (b)). In this case one 
expects that the Coulomb interaction between electrons 
located in the active states causes very small changes in 
the transient currents. On the other hand, the Coulomb 
repulsion generated by the inactive states which are fully 



occupied could be important, leading at least to a Hartree 
shift of the active levels. However, the relaxation pro- 
cesses and the modulation of the output current should 
be qualitatively similar. 



IV. CONCLUSIONS 

We have analyzed the transient response of a two- 
dimensional nanosystem to a sequence of periodic rect- 
angular pulses which modulate the contact to the source 
lead. By solving the generalized non-Markovian master 
equation for the reduced density matrix we have been 
able to discuss the dependence of the input /output char- 
acteristics on the pulse length for two specific systems: 
a rather large quantum dot and a narrow quantum wire. 
We have considered a pump-and-probe setup in which the 
contact to the left lead opens and closes periodically while 
the right lead is always connected to the system. When 
the contact is switched off the drain current reflects the 
relaxation processes in the sample. We have discussed 
these processes by analyzing the single-particle currents, 
the diagonal, and off-diagonal elements of the reduced 
density matrix. At a low temperature the phonon effects 
could be neglected and the main relaxation processes are 
back-and-forth tunnclings to and from the leads. 

In both cases (the dot and the wire) we have found that 
the pulse length can be adjusted such that the shape of 
the pulse can be reproduced by the output signal. By 
increasing the chemical potential of the source lead the 
current profile in the output lead develops additional os- 
cillations related to the relaxation of the higher energy 
states included in the bias window. In the case of a nar- 
row quantum wire a delay of the output current with 
respect to the input signal has been obtained. 

This study was partly motivated by the recent exper- 
iments of Naser et al^ and Lai et alA Although those 
experiments were done with larger and longer pulses our 
results qualitatively agree with the reported features of 
the transient response. In particular the time-dependent 
occupation number clearly show a charging/relaxation 
behavior as in the work of Lai et alA 
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